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Abstract 

When dealing with certain kind of complex phenomena the theoretician may face some difficul- 
ties - typically a failure to have access to information for properly characterize the system - for 
applying the full power of the standard approach to the well established, physically and logically 
sound, Boltzmann-Gibbs statistics. To circumvent such difficulties, in order to make predictions 
on properties of the system and looking for an understanding of the physics involved (for exam- 
ple in analyzing the technological characteristics of fractal-structured devices) can be introduced 
large families of auxiliary statistics. We present here a systematization of these different styles 
in what can be termed as Unconventional Statistical Mechanics, accompanied of an analysis of 
the construction and a clarification of its scope and interpretation. As illustrations are derived 
heterotypical Bose-Einstein, Fermi-Dirac and Maxwell-Boltzmann distributions, and several appli- 
cations including studies of experimental works are briefly described. 

PACS numbers: 05.70.Ln, 82.20.Mj, 82.20.Db 

Keywords: Nonequilibrium Statistical Mechanics; Nonextensive Thermo-Statistics; Renyi Statistics; Uncon- 
ventional Statistics; Informational Entropies, Escort Probability, Fractal-Structured Systems 



1 



I. INTRODUCTION 



Over twenty years ago Montroll and Shlesinger wrote that in the world of investigation of 
complex phenomena that requires statistical modeling and interpretation, several competing 
styles have been emerging, each with its own champions Q]. Lately a large amount of efforts - 
with a flood of papers related to one particular case - has been given to the topic. What is at 
play consists in that in the study of certain physicochemical systems we may face difficulties 
when handling situations involving fractal-like structures, correlations (spatial and temporal) 
with some type of scaling, turbulent and chaotic motion, small size (nanometer scale) system 
with eventually a low number of degrees of freedom and complicate boundary conditions, 
generalized hydrodynamics, and so on. The interest on them has been recently enhanced 
as a consequence that such situations are present in electronic and opto-electronic devices 
of the nowadays point-first technologies and also in technological areas involving the use 
of disordered systems, polymeric solutions, conducting glasses, the case of microbatteries, 
and others. The question is that these difficulties consist, as a rule, in that the theoretician 
cannot properly satisfy Fisher's criteria of efficiency and sufficiency 52| in the conventional, 
well established, physically and logically sound Boltzmann-Gibbs statistics, meaning an 
impairment to correctly include the presence of large fluctuations (and eventually higher- 
order variances) and to account for the relevant and proper characteristics of the system, 
respectively. Then, just out of necessity and convenience and to be able to make predictions 
- and in the way providing an understanding, even partial, of the physics of the system of 
interest, for example in analyzing technological characteristic of physico-chemical systems 
as illustrated in Sect. Ill - one may resort to the use of statistics other than the usual form 
of the Boltzmann-Gibbs (BG) one, which are not at all extensions of the latter but, as said, 
introduce an auxiliary detour for analyzing the problem in hands. 

Among existing approaches we can mention what can be called Generalized Statistical 
Mechanics as used by P. T. Landsberg showing that functional properties of the (informa- 
tional) entrop,es (see next secfon) give or,g,„ to efferent types of thermostatisties. and r.se 
the question of how to select a "proper" one, that is, are some better than others? |3| ; for 
decades has been in use Levy Statistics, introducing a modified non-Gaussian distribution 
which has been applied to a variety of problems (see for example Ref. 0]); the approach of 
W. Ebeling who has addressed the question of the treatment of a class of systems, in nature 
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and society, which are determined by their total history; it can be referred-to as Ebeling 
Statistics p; Super statistics developed by C. Beck and E. G. D. Cohen for nonequilibrium 
systems with complex dynamics in stationary states with large fluctuations on long-time 
scales !&]; Non-Extensive Statistics based on Havrda-Charvat statistics 3| applied to a num- 
ber of cases are described in the Conference Proceedings of Ref. j8|: It consists in the 
so-called Tsallis statistics j^; Renyi statistics - devised by the renowned theoretical statis- 
;itian A. Renyi llyl - has been introduced in the scientific literature, as noticed in Ref. 

~l n 

with, for example, P. Grassberger and I. Procaccia 12] using it as a valuable method 
for characterizing experimental chaotic signals, and recently P. Jizba and T. Arimitzu 
have presented an extensive analysis of it in a paper called "The world according to Renyi" ; 
Sharma-Mittal Statistics or better to say a variation of it (called Kappa or Deforma- 
tional statistics) was used by V. M. Vasyliunas in problems of astrophysical plasma [15] and 
by G. Kaniadakis in the case of relativistic systems and so on. 

It has been stated that some of these statistics are nonextensive as compared with BG 
statistics (e.g. Ref. j^) but this is not quite correct once in the strict equilibrium of 
adiabatically closed systems, described by the microcanonical ensemble, the BG entropy is 
nonextensive only becoming practically extensive in the thermodynamic limit (e.g. Ref. 17[). 
The latter is one case of idealization in science hsf, useful for particular situations however 

n 

not applicable in general: notice for example the case of small systems [19|| present, e.g., in 
nano-science and technology which are of a quite large interest and in ample development 
nowadays. 

We recall that in Statistical Mechanics the probability distribution (statistical operator), 
usually derived from heuristic arguments, can also be derived from a variational method, once 



is made connection with Information Theory 
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22| . It consists into making extremal - 



a maximum -, subject to certain constraints, a functional (superoperator) of the probability 
distribution (statistical operator). Such quantity, first introduced in Shannon's Theory of 
Communication j22i, can be referred to as measure of uncertainty of information. It has also 
been called statistical measure and entropy, with the understanding that it is informational 
entropy. It is worth to emphasize - in view of some confusion that has recently pervaded 
the scientific literature - that the different possible informational entropies are not to be 
interpreted as the thermodynamic entropy of the physical system. R. T. Cox has noticed 
that the meaning of such entropies is not the same in all respects as that of anything which 
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e to give a simple verbal 



2J|. E. T. Jaynes has also 



has a familiar name in common use, and it is therefore impossib 
description of it, which is, at the same time, an accurate definition 
commented that it is an unfortunate terminology, and a major occupational disease in that 
there exists a persistent failure to distinguished between the informational entropy, which is 
a property of any probability distribution, and the experimental entropy of thermodynamics, 
which is instead a property of the thermodynamic state: Many research papers are flawed 
fatally by the authors' failure to distinguish between these entirely different things, and 
in consequence providing nonsense results 25|. Along such erroneous line, recently it has 
been considered the utterly wrong proposition that a particular informational entropy - 
among the infinitely-many that can be defined - comes to supersede the Boltzmann-Gibbs 
one as the entropy of physical systems 0, : Such "entropy" has a form of the so-called 
structural Havrda-Charvat one j7|, which is, we restate, a generating functional for deriving 
so-called heterotypical probability distributions js?! to contour, as noticed, the difficulties of 
application of Boltzmann-Gibbs statistics when we face limitation (ours, not of EG statistics) 
in satisfying Fisher's criteria in statistics j^]. 

As noticed in the Abstract we present a systematization of these alternative statis- 
tics, dubbed Unconventional Statistical Mechanics (USM), within the variational formalism 
MaxEnt-NESEF (for Maximization of Informational Entropy in the Non-Equilibrium Sta- 
tistical Ensemble Formalism j^^Sl) with the use of non-conventional statistical measures 
(informational entropies) to derive heterotypical nonequilibrium probability distributions. 
This is discussed in the next Section, and after that we illustrate the matter with the deriva- 
tion of heterotypical Fermi-Dirac, Bose-Einstein, and Maxwell-Boltzmann single-particle 
distribution functions, and with a brief description of the application of USM in several 



situations including analysis of experimental data (see also Refs. 
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lad). 



II. UNCONVENTIONAL STATISTICAL MECHANICS 



In Statistical 
Introduction - 



Mechanics the variational approach MaxEnt-NESEF - mentioned in the 
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provides a powerful, practical and soundly based procedure, 
of a quite broad scope, for building a nonequilibrium ensemble formalism. It is encompassed 
in what is sometimes referred-to as Informational- Entropy Optimization Principles (see for 
example Ref. ^^), or, to be more precise, we would say constrained optimization, that 
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is, restricted by constraints consisting in the available information. Such optimization is 
performed through calculus of variation with Lagrange's method for finding the constrained 
extremum being the preferred one. 

In the conventional approach, let it be in the cases of equilibrium js^ls^' nonequilib- 
riurn (in the linear or Onsagerian regime 3^ , or for arbitrarily far-from-equilibrium systems 



mm), 

one proceeds to maximize Shannon- Jaynes measure of information - or infor- 



mational entropy which in equilibrium is Gibbs entropy giving rise to Boltzmann entropy in 
the microcanonical ensemble in equilibrium -, namely 



SsAt) = -Tr{pit)lnpit)}, (1) 

under the given constraints. The calculus leads to statistical operators of an exponential 
form, with the exponent depending on the set of basic dynamical variables present in the 
constraints {{Pj{r)}, see below), and the set of Lagrange multipliers {{Fj{r,t)}, see below), 
the latter being the nonequilibrium (or equilibrium according to the case) thermodynamic 
variables. In fact we must proceed to the maximization of the SJ-informational entropy with 
the constraints 



Tr{p(t)} = l 
Q,(r,t)=Tr|p,(r)p(t) 



(2) 
(3) 

with J = 1,2, . . .. Equation Q is the normalization condition and Eq. consists in the 
set of average values over the nonequilibrium ensemble (characterized by the sought after 
statistical operator p(t)) of the set of microdynamical variables {Pj{r)} chosen for providing 
;he constraints in the variational method. The resulting statistical operator is given by 
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Pe(t)=exp|-5(t,0) + dt'e<''~''>^S{t',t'-t)Y 



(4) 



where 



S{t',t' -t) =exp{- — {t'-t)HjS{t',0)exp^—{t'-t)H}, 



(5) 
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and 



(6) 



is the so-called informational-statistical-entropy operator [37[, and p{t,0) is an auxiliary 
(sometime called "instantaneously frozen quasiequilibrium" ) statistical operator having an 
exponential form resembling a canonical-like distribution 29] . In these expressions, H is the 
system Hamiltonian, {Fj{r,t)} the set of Lagrange multipliers associated to the set of basic 
microdynamical variables {Pj{r)}, and (j){t) ensures the normalization condition of Eq. (jSJ 
and can be considered as being the logarithm of a nonequilibrium partition function, say, 
(f)[t) = lnZ{t). Moreover, the term containing the positive infinitesimal e - which goes to 
zero after the calculation of averages has been performed - results from introducing Abel's 

We 



kernel (in the convergence of integral transforms) in the integral on time |22l . |2 
recall that it introduces in the theory the concept of Bogoliubov's quasiaverages [38| leading 
to irreversible evolution from an initial condition, what it does by selecting the retarded 
solutions of the Liouville equation that p{t) satisfies: the advanced solutions are discarded 
in a quite similar way as done by Gell-Mann and Goldberger in the case of Schrodinger 

satisfies, as it should, 
4l| - and irreversibility 



equation in scattering theory [39[ . The statistical operator of Eq. 
the question of historicity - in Kirkwood's and Mori's sense |4(, 



is introduced here in the so-called interventionistic picture in logic, resorting to Krylov's 

using a Poissonian distribution (Abel's kernel) to account for 



"jolting" approach |4 



22 



33|. 



the nonisolation of any physical system in the real world 

In that way it is built a general ensemble formalism in Boltzmann-Gihhs Statistical Me- 
chanics for any kind of thermodynamic state of the system, and fundamental for the study 
of the physics of condensed matter: it allows to describe successfully any possible physical 
situation in many-body systems. Of course, such success is possible if we are able to prop- 
erly handle it, what means, as already mentioned, that in any problem we need to satisfy 
Fisher's criteria ^ for it to give satisfactory predictions. As noticed, we may face difficul- 
ties for some kind of situations where the use of BG statistics is simply impaired because of 
either becomes not possible to handle the required information relevant to the problem in 
hands, or, more important and fundamental, is involved a failure on the part of the theoreti- 
cian to have a correct access to such information. Typical cases are that of small systems 
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(few particles) when we may not be satisfying Fisher's criterion of efficiency (fluctuations 
and higher-order variances cannot be ignored), and the case of complex systems with some 
type or other of fractal-like structures or with long-range space correlations or particular 
long-time correlations, when we have deficiencies in the proper knowledge of the character- 
ization of the states of the system (failure to satisfy Fisher's criterion of sufficiency). We 
reemphasize that what is present is a practical difficulty (a limitation on the part of the 
theoretician) in the otherwise complete, physically and logically sound BG-statistics, and 
this is illustrated in Section III. 

When facing such difficulties, out of necessity (obstacles in accounting for the relevant 
information) or, mainly, the force of circumstances (lack of the proper characterization of 
the system) a way to circumvent them has consisted into resorting to auxiliary (or uncon- 
ventional) statistics (see the first sentence in the Introduction). 

This auxiliary Unconventional Statistical Mechanics consists of two steps: 

1. The choice of a heterotypical probability distribution, obtained through application of 
MaxEnt to different measures of information, that is, other than the Shannon- Jaynes 
one that leads to the conventional BG-statistics; 

2. The use of the escort probability in terms of the chosen heterotypical probability of 
item 1 above, in the calculation of averages. 

Let us consider both concepts. The escort probability of order 7 of a given probability 
distribution p is defined as 

V,{p}=pyTr{p^}, (7) 

where 7 is a positive real number, which seems to have been introduced by A. Renyi but in 
terms of Renyi's heterotypical probability distribution (see Ch. 10 in the book of Ref. Q), 
being generalized, and the name given, by C. Beck and F. Schlogl (see Ch. 5 in the book 
of Ref. 13 )• ^ol^ is normal definition of average value the presence of 

second and higher-order variances (in that way trying to improve upon the possible failure 
of distribution p to satisfy Fisher's criterion of efficiency): For the average value of an 
observable A, if we write 7 = 1 + e and make an expansion in powers of e we find that 
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i) = Tr[Av,{p}} 



e 

+ 2 
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2(A)(Sy -2/as\(S 



(8) 



where 



(...)=Tr{...p}, (9) 

stands for the normal average. For illustration let us take for A the Hamiltonian H and for 
p a canonical distribution p = Z^^ exp {—PH}, and then using Eq. (jH)) up to second order 
in e = 7 — 1 it follows that the energy is given by 



E 



H 



H)+ePA2E + -P'AsE, 



where 



(10) 



A^E 
AsE 



H -IH 



H -{H 



H 



3 H) H 



2(H 



(11) 
(12) 



with A2E and A^E being the second and the third order variances of the energy. The 
departure of 7 from the value 1 gives an indication of the influence of the variances in the 
process of prediction of average values. 

On the other hand, for obtaining the heterotypical distributions one needs to introduce 
alternative measures of information (i.e. informational entropies). A large family is the one 
provided by I. Csiszer ([4^ and see also the book in Ref. 27 1) and, among the infinitely-many 
possibilities, we list in Table I five of them (M^~^ is the constant probability in the uniform 
distribution). Kiilback-Leibler measure 4^, which is parameter independent, corresponds 
to Shannon- Jaynes informational entropy and produces the usual results in BG-statistics. 
The others are dependent on what has been dubbed as infoentropic index(es): The one 
of the Havrda-Charvat Q] has been used in Physics and other disciplines with the name 
of Tsallis entropy (calling q the infoentropic index); Sharma-Mittal informational en- 
tropy [1^ depends on two infoentropic indexes (taking index /3 = 1 it goes over the one of 
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Havrda-Charvat, and for a particular relation between both indexes results Kappa-statistics 



of Refs. 



which is a kind of weighted superposition of two Havrda-Charvat mea- 



sures. Renyi measure 10 [has been used in a number of problems in several disciplines, 
including Physics jlll llj . isl ] (it can be noticed that Renyi and Havrda-Charvat measures 
produce, via MaxEnt, the same statistical operator, differing only in the interpretation of 
the Lagrange multipliers). Kapur measure 47[, depending on two infoentropic indexes is a 
superposition of two Renyi measures. 

Two important points need be stressed: On the one hand, Shannon-Kiilback-Leibler mea- 
sure is derived from a set of fundamental axioms, and Havrda-Charvat and Renyi measures 
follow from the modification of one of such axioms 0, 10, 46 1. On the other hand, in the 
process is introduced an open parameter (the infoentropic index a, recently indicated by q), 
to be determined by fitting of the index-dependent prediction with observation. 

The role of these heterotypical probability distributions is to introduce a selective 
weighting of, for example, occupation functions of quantum states, amplitudes of movement 



in hydrodynamics, etc. This is illustrated in next Section. 



Table I: Informational-Statistical Entropies 



Conventional ISE 


Boltzmann-Gibbs-Shannon-Jaynes ISE 
(from Kiilback-Leibler measure) 


~Tr{p In p] 


Unconventional (entropic-indexes-dependent) ISEs 


From Havrda-Charvat measure 


-^Tr{p" - p] 
a > and a ^ 1 


From Sharma-Mittal measure 


a>l, l3<loia<l, 13>1 


From Renyi measure 


a > and a ^ \ 


From Kapur Measure 


-^[lnTr{p''] - InTripf^}] 
a > 0, /3 > and a ^ {3 
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III. ILLUSTRATIVE EXAMPLES 



In this Section we illustrate the use of the Unconventional Statistical Mechanics described 
in the previous ones. We present comments on a few cases and is given the reference 
where details can be consulted. One example is the case of "anomalous" luminescence 
in semiconductor heterostructures, which requires the introduction of heterotypical Fermi- 
Dirac distributions. We begin with this case. 

A. Photoluminescence in Semiconductor Heterostructures 

There exists nowadays a large interest on the question of optical properties of quantum 
wells in semiconductor heterostructures, which have been extensively investigated in the last 
decades as they have large relevance for the high performance of electronic and optoelec- 
tronic devices (see for example Ref. j^)- ^^^^ with these kind of systems, because of 
the constrained geometry that they present (where phenomena develop in nanometer scales) 
the researcher has to face difficulties with the theoretical analysis: A most relevant ques- 
tion to be dealt with is the one related to the interface roughness, usually having a kind 
of fractal-like structure which leads to energies and wavefunction depending on boundary 
conditions which need account for spatial correlations. As a consequence the different phys- 
ical properties of these systems appear as, say, "anomalous" when the results are compared 
with those that are observed in bulk materials, particularly, the case of photoluminescence 
which we briefly describe here. In the study of photoluminescence in nanometric quantum 
wells the conventional treatment via the well established Boltzmann-Gibbs formalism has 
its application impaired because of the spatial correlations resulting from the non smooth 
confinement, relevant in the characterization of the system, on which one does not have 
access to (obviously the interface roughness varies from sample to sample and one does not 
have any easy possibility to determine the topography of the interface). This is then the 
reason why the criterion of sufficiency is not satisfied in this case. 

Let us consider a system of carriers (electrons and holes) produced, in the quantum well 
of an heterostructure, by a ultrafast laser pulse. They are out of equilibrium and their 
nonequilibrium macroscopic state can be described in terms of a nonequilibrium statistical 
thermodynamics with the nonequilibrium thermodynamic state characterized by the 
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time evolving quantities energy and density (see for example Ref. Electrons and holes 

do recombine producing a luminescence spectrum which, we recall, is theoretically expressed 

as 

OC Yl /nkj^) fn'^A^)^ - - ^nkj, (13) 

n,n',kx 

where and are the populations of electrons and holes, HQ = huj — Eg, with uj being the 
frequency of the emitted photon and Eg the energy gap, and e^^]^| are the electron (hole) 
individual energy levels in the quantum well (index n for the discrete levels and for the 
free movement in the x — y plane), and the delta function accounts for energy conservation. 

The textbook expression for the energy levels corresponding to the use of perfectly smooth 
bidimentional boundaries is given by 

where Lqw is the quantum-well width and '^e(/i) effective mass and the populations 

/ take a form that resembles instantaneous in time Fermi-Dirac distributions: see Ch. 6 in 
the book of Ref. j^] and Ref. jj49| . Using this expression in the calculation of / of Eq. ()13|) 



and comparing it with the experimental results |5Cll. Isi | one finds a disagreement, and then 
it is used the name of "anomalous" luminescence for these experimental results. This is a 
consequence that we are using an improper description of the carriers' energy levels - we 
are not satisfying the criterion of sufficiency (as discussed in previous Sections) resulting 
of ignoring the roughness of the boundaries (with self-afiine fractal structure j53]) which 
needs be taken into account in these nanometric-scale geometries, and then the boundary 
conditions to be placed on the wavefunctions are space dependent. Hence complicated space 
correlations are to be introduced, but to which we do not have access (information), as 
already noticed. Hence, this limitation on the part of the researcher breaks the sufficiency 
criterion, and application of the Boltzmann-Gibbs-Shannon-Jaynes construction is impaired, 
and one can try to circumvent the difficulty introducing unconventional statistics based on 
parameter-dependent structural informational entropies in order to analyze the properties 
of the system. 
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1. Distribution Fuctions of Fermi ons and Boson in USM 



According to Eq. p3p what we need is to express the carrier populations in a heterotypical 
statistics, and we resort to the Renyi one (see for example Ref. 



Using 



lobulations ii 
PJ and laO 
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MaxEnt-NESEF in terms of Renyi statistical entropy it can be derived the corresponding 
auxiliary ("instantaneously frozen") statistical operator p(t, 0) [cf. Eq. ® - the case of 
the conventional one and we recall that the statistical operator is given in term of this 
auxiliary one in Eq. (0))]. After some straightforward mathematical manipulations it follows 
that it can be written in a convenient form for performing calculations, namely. 



Pa{t,0) 



Vait) 



[a 



1 

' a-l 



(15) 



where 



fjait) = Tr 



1 + {a 



[a 



(16) 



(17) 



Equation stands for a modified form of the quantity that ensures the normaliza- 
tion condition, and Eq. (fT7j) for modified Lagrange multipliers, where Fma are the original 
Lagrange multipliers, and the average values of the microdynamical variables Pm- More- 
over, we recall that in Renyi statistics the associated escort probability is given by 



(18) 



that is, the order of the escort probability is the same as the index in Renyi heterotypical 
distribution (see Ch. 10 in the book of Ref. Ql)- In particular to derive the distribution 
functions for fermions and for bosons using USM in terms of Renyi statistical approach one 
chooses as basic dynamical variables, i.e. the Pj in Eq. ()15|). the set of occupation number 
operators 



{nk} = {c^Ck}, 



(19) 



12 



where c{c^) are the usual annihilation (creation) operators in states |k), satisfying the cor- 
responding commutation and anticommutation rules of, respectively, bosons and fermions 
(the spin index is ignored). Their average values are the infoentropic-index ct-dependent 
distribution functions 



/k(t) = Tr |4ck2^,,{Pae (t)}} = Tr |4ckPa{pa(t, 0)}} 



(20) 



which is a consequence that for the basic variables, and only for the basic variables, the 

average value coincides with the one 

auxiliary statistical operator is then [cf. Eq. (|15|)] 



Pa{t,0) 



1 - (a 



1) J]^ka(t)4ck 



(21) 



with [cf. Eq. (p7|l ] 



n -1 



l + (a-l)5^Fk',(t)/k'(t) • (22) 

n 

The populations of Eq. (|^. according to the calculation described in Ref. j30| take the 
form 



/k(t) = /k(t)+Ck(t), 



(23) 



where 



1 + (« - l)Fka(t) 



(24) 



± 1 



where upper plus sign stands for fermions, and the lower minus sign for bosons, and 



Ck(t) = a{l - a)(l - h{t))J2ha{t)Fi,,^{t)Tr {clc^clc^'VMt,^)}} + ■■■, (25) 

k 

involving two, three, etc. particle correlations, which in general are minor corrections to the 
first, and main, contribution, the one given by Eq. ()24j) . 
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In the limit of a going to 1, which apphes when the criteria of efficiency and sufficiency 
are satisfied, Renyi statistical entropy acquires the form of the Boltzmann-Gibbs-Shannon- 
Jaynes one, C becomes null, -fka(i^) becomes Fk(^), and then 

AW = ^Jyibl- 

(In equilibrium -Fk(^) (^k — l^l^B^ and there follows the traditional Fermi-Dirac and 
Bose- Einstein distributions) . 

We can see that the distribution of Eq. ()23|) is composed of a term / corresponding to 
the individual particle in state |k), plus the contribution C containing correlations (of order 
two, three, etc.) among the individual particles. This type of calculation but for systems 
in equilibrium and not using the average value defined in Eq. (jH)), in terms of the escort 
probability, was reported in Ref. js^, and then are not satisfactory. 

Let us now give some attention to the Lagrange multipliers -Fka(^)- The most general 
statistical operator for nonequilibrium systems can be expressed in the form of a generalized 
nonequilibrium grand-canonical statistical operator for a system of individual quasiparticles, 
where the basic variables are independent linear combinations of the single-quasiparticle 



occupation number operators fc 



and their fiuxes of all order 
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". Eq. (HD)], consisting of the energy and particle densities 
29|. This description follows from the choice 



Kaif) = /3«(t)[ek - fia{t)] - Fha{t) ■ eku(k) - F„,(t) • u(k) 

- E [^'iW ® ^i^"'^'(k) + Fj[^{t) ® t.M(k)] , (27) 



r>2 



where has been introduced the scalar quantities P(t) and fiait), Fha(t) and F„Q,(t) are vectors. 



and and Fnd r-th rank tensors. Moreover, 



MM(k) = [u(k) . . . (r - times) . . . u(k)] , (28) 

is the tensorial product of r-times the characteristic velocity u(k) = h'^V^ek, where ek 
is the energy dispersion relation of a single-particle, and then u(k) is the group velocity 
in state |k). Dot stands as usual for scalar product of vectors, and ® for fully contracted 
product of tensors. 
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To better illustrate the matter, we introduce a simplified description, or better to say a 
quite truncated description, proceeding to neglect in Eq. (jTTj) all the contributions arising 
out of the fluxes, i.e. we put F = and -F''"' = 0, retaining only the first term on the 
right-hand side. Therefore, we do have that 



hit) = , (29) 



1 + {a - l)f3ait)[ek - fiait)] 



where [cf. Eq. 



Pa{t) = f3ait)/{l -{a- l)(3^{t)[E{t) - fiait)N{t)]} , (30) 

and Pfl = Pfi, then fl = having the role of a quasichemical potential. 
In this Eq. (IHUj) E(t) is the energy 

E(t)^^ek/k(t), (31) 

k 

and the number of particles 

iV(t)^5^/k(t), (32) 

k 

where the correlations present in C in Eq. (PHj) have been neglected. Moreover, in many cases 
we can use an approximate expression for the populations, that is, in the one of Eq. 
we admit that ±1 can be neglected in comparison with the other term. This is considered 
as taking a statistical nondegenerate limit, once, if we put a going to 1 (what, we again 
stress, strictly corresponds to the situation when the principle of sufficiency is satisfied), 
the population takes the form of a Maxwell-Boltzmann distribution. In this condition the 
expression for the population can be written as 

Mt) = A„(t)[l + (« - l)i3„(t)ek]-^ , (33) 

where 

A,(t) = [l-{a- l)(3a{t)fia{t)r^^ , (34) 



Bo.{t) = Pa{t)/[l -{a- l)/3„(t)/i«(t)] 
15 



(35) 



Consider a parabolic dispersion relation, that is, ek = h'^k'^ /2m* . Using Eq. ()33p in Eqs. 
(p?T|l and (jn21), we arrive at the result that 



n{t) = ^ = ATit)^hda), (36) 



.,/3/2(a), ^ 

e{t) = —— = n{t)——-kB%{t) , (37) 

V h/2[(y) 



with the integrals lu^a) related to Beta functions 0|, and we have introduced the definition 

B-\t) = kB%{t) , (38) 

where T plays the role of a pseudotemperature and where in Eq. ()36|) is a characteristic 
length given by A^(t) = fi^ /m*kBTa{t) (that is, de Broglie wave length for a particle of mass 
m* and energy kBTaif)). 

We can see that the above Eqs. and (fHTj) define the Lagrange multipliers, /3Q,(t) and 
jiait) - present in Aaif), Xa{t), and Bait) - in terms of the basic variables energy and number 
of particles. Moreover, using Eq. (jHUj) we can obtain an expression for the quasi-chemical 
potential in terms of quasitemperature and density, namely 

1 - (« - l)/3„(t)/i,(t) = [47i'\l{t)/h/2{a)]''^ [n{t)f-^ . (39) 

Also, it can be noticed that for a = 1 (the case when the condition of sufficiency is 
satisfied) one recovers the equivalent of the results of conventional nonequilibrium statistical 



mechanics 
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36|, which are 



e{t) = ^n{t)kBT*{t), (40) 

where we have introduced the so-called quasitemperature defined by fc^T* = 

i3~l^(t), this equation standing for a kind of equipartition of energy at time t, and 

/i(t) = -kBT*{t) ln[T*it)/etr{t)] , (41) 

where 9tr{t) = ^^n^/^(t)/2m* is the characteristic temperature (here in nonequilibrium con- 
ditions and at time t) for translational motion. This suggests to define a so-called "kinetic 



16 



temperature" Qxit) |57[ by equating e(t) to {3/2)n(t)kBQ xit) , given, after Eq. (|H7j) is used, 
by 



0^(t)=T,(t)/(5-3a), (42) 

where we can see that a must be smaller than 5/3, as shown in Ref. j^^] where connection 
of theory with experiment is presented, together with other illustrations and discussions. 

Let us consider how does the a-dependent distribution of Eq. compares with the 

usual Fermi-Dirac and Bose-Einstein distributions. For illustration we consider the nonde- 
generate limit of Eq. (j33p . common to both, where parameter B is related to the kinetic 
temperature Qk by Eqs. and (ji^ . Taking for 7^ of Eq- ffHH]) the unique value of 

3007^, Figs. [U and El show the population of Eq. (jH!^ corresponding to several values of the 
infoentropic index a. It can be noticed that the formalism introduces the characteristic of 
a different weighting of the values of the standard distribution {a ~ 1), such that: (1) for 
a < 1 the population of the modes at low energies are increased at the expense of those of 
higher energies (e > 7 x lO'^^eV), while (2) for a > 1 we can see the opposite behavior. 



2. The "Anomalous" Photoluminescence Spectrum 

Returning to the question of "anomalous" luminescence in nanometric quantum wells in 
semiconductor heterostructures, the carrier's populations, to be used in Eq. ((131), of ^be 
form of Eq. , namely 

}^±l}", (43) 

which depend on e^^]^|, that is, the ideal single-carrier energy level of Eq. ()14|). Using these 
populations in the nondegenerate limit [cf. Eq. ^^], the luminescence spectrum of Eq. 
f|T3|l is given by 




I{uj\t) oc 



m 



ml 



m 



ml 



1 + (a - mffM^n + (a - lYBim^^it) 



m„ 



(44) 
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where f3efu = ^Bl{t) + ^S^(t), m'^ = [m:]-i + [ml]-\ and M = ml + ml. The 



second contribution in the last term in Eq. (jl^ is much smaller than the first, as verified a 
posteriori, and then the spectrum is approximately described by 

I{uj\t) (X [1 + (« - l)^ffSt){huj - Eg)]^- . (45) 
The experiments reported in Ref. are time integrated, that is, the spectrum is given by 

X{uj) = —j^ dt'Iiu;\t') , (46) 

where At is the resolution time of the spectrometer. Using Eq. ()45p in Eq. ()46|) . and in the 
spirit of the mean- value theorem of calculus we write 

I{uj) oc [1 + (a - l)^effA^uj - Eg)] ^ , (47) 

introducing the mean value Peff^ (as an open parameter), which we rewrite as [Peffa]^^ = 
kB^ai defining an average, over the resolution time At, effective temperature of the nonequi- 
librium carriers, that is, a measure of their average kinetic energy (see Ref. js?!). 

In Fig. ini is shown the fitting of the experimental data with the theoretical curve as 
obtained from Eq. (|47p . It contains the results referring to four samples having different 
values of the quantum well width. The information-entropic index a depends, as expected, 
on the dimensions of the system: as the width of the quantum well increases the values of a 
keep increasing and tending to 1 (see Fig. |3I). This is a clear consequence that the fractal- 
like granulation of the boundary surface becomes less and less relevant for influencing the 
outcome of the phenomenon, as the width of the quantum well falls outside the nanometer 
scale, and is approached the situation of a normal bulk sample. On the other hand the kinetic 
temperature of the carriers is smaller with increasing quantum well width, as also expected 
once the relaxation processes, mainly as a result of the interaction with the phonon system, 
become more effective and the cooling down of the hot carriers proceeds more rapidly. 

Moreover, it can be empirically derive what we term as a law of " path to sufficiency", 
namely, 

^ (48) 
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where, by best fitting, Li ~ 139 ± 17, L2 — 204 ± 24, all values given in Angstrom. We do 
have here that as L largely increases, the entropic index tends to 1, when one recovers the 
expressions for the populations in the conventional situation, but as L decreases a tends 
to a finite value L1/L2, in this case ^ 0.7 ± 0.06. This indicates that the insufficiency of 
description when using Eq. (fT^ in the calculations (the ideal energy levels) becomes less 
and less relevant as the size of the system increases as already commented. 

This illustration of the theory clearly evidences the already stated fact that the infoen- 
tropic index a is not a universal one for a given system, but it depends on the knowledge of 
the correct dynamics (the region of energy-momentum space that is involved in the exper- 
iment being analyzed), the geometry and size including the characteristics and influence of 
the boundary conditions (e.g. the fractality in the present case), the macroscopic (thermody- 
namic) state of the system in equilibrium or nonequilibrium conditions, and the experimental 
protocol. 

The case we presented consisted of experiments in time-integrated optical spectroscopy. 
The phenomenon of "anomalous" luminescence in nanometric quantum wells in semicon- 
ductor heterostructures is also present in the case of time-resolved experiments (nanosecond 
time resolution where the infoentropic index and the (nonequilibriurn) kinetic temperature 
change in time accompanying the irreversible evolution of the system pS'). The use of USM 
in a completely analogous way as done above (note that Eq. is valid for a time-resolved 
situation) allows to determine the evolution in time of the kinetic temperature 6Q,(t), and 
the infoentropic index a(t), which is then changing in time as it accompanies the evolution 
of the irreversible processes in the nonequilibrium thermodynamic state of the carriers. 

This case involves, as noticed, the difficulty of insufficiency in the characterization of 
the energy levels of the carriers in the quantum well, i.e. it is of a microscopic mechanical 
character. As shown before, use of the escort probability introduces ad hoc information 
trying to account for missing relevant correlations in the problem, and Renyi (heterotypical) 
statistical operator produces modifications in the carriers' distributions trying to account 
for insufficiency of knowledge about the proper quantum mechanical levels. 



19 



B. "Anomalous" Diffusion 

On the other side, insufficiency of description at a macroscopic level is present in ques- 
tions involving hydrodynamics. So called "anomalous" situations are also presented here, 
e.g. "anomalous" diffusion (non-Fickian diffusion) 5^ . The question now is which is the 
source of the lack of sufficiency in the well established Boltzmann-Gibbs formalism, which 
forces us to resort to the unconventional approaches. The answer resides in that is being 
used a quite incomplete hydrodynamic approach in situations which require a more detailed 
treatment. In its more general approach, the description of hydrodynamical (including 
rheological) motion should consist of an extended description, which can be referred-to as 
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62| . Ignoring fluctuations 



Non-Linear Higher-Order Hydrodynamics with fluctuations 
(relevant for example in turbulent motion) one needs to introduce the densities of energy, 
h{r,t), and particles, n(r, t), and their fluxes of all order, namely Ip(r, t), for the first (vec- 
torial) fluxes, Ip\r,t) for the higher-order (r-rank tensor) fluxes, r > 2, and p = h or n 
for energy and particle (material) motion respectively. The motion is then determined by a 



complicated set of equations of evolutions of the type jZ 
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62, Q 



§-/;Kr,t) + V- (r , t) = (r, t) , (49) 

where r = for the density, r = 1 for the first (vectorial) flux, or current, and r > 2 for 
the all higher-order fluxes, J7n are collision integrals, and V- is the operator indicating 
to take the divergence of the tensor. Solving this set of coupled equations of evolution is a 
formidable, almost unmanageable, task. As a rule one uses, depending on each experimental 
situation, a truncation on the set of equations (i.e. they are considered from r = up to a 
certain value, say n, of the order r) Hence, if one is restricted to introduce a low order 
truncation it is faced a failure of the criterion of sufficiency when using the conventional, and 
universal, approach. Consequently, a way to circumvent the difficulty is, as shown before, to 
make calculations in term of unconventional statistical mechanics, e.g. resorting to Renyi's 
approach. 

The necessity to go over higher orders comes, for example, because of the presence of 
a fractal-like structure in the system. We have dealt in detail the question of polymeric 
solutions 6^, where the question is presented and discussed in detail being shown how 
there follows in USM "anomalous" Maxwell-Cattaneo and diffusion equations. In this case. 
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macromolecules under flow, one faces the difliculty in the description resulting from the 
self-similarity in a type of average fractal structure that this shows (what we have called 
"Jackson Pollock Effect" , in view of the analogy with his paintings with the dripping method, 

given elsewhere ^6^, the 
"anomalous" diffusion equation, derived in Renyi statistics, has the form 



^n(r, t) - D^{r, t)V'n^- (r, t) = -r,(r, t) V ■ V ■ (n(r, t) [v(r, t)v(r, t)]) + ... , (50) 

where Da{r,t) is a transport coefficient whose dependence on position and time is a con- 
sequence of its dependence on the nonequilibrium thermodynamic state of the system, and 
so is the case of the momentum relaxation time Ta{r,t); moreover, the first term on the 
right - depending on the velocity field - is the double divergence of the so-called convective 
pressure tensor ([vv] stands for tensorial product of twice the velocity vector rendering a 
rank-2 tensor), and we have omitted to write down additional terms involving gradients of 
the transport coefficient and the momentum relaxation time. Finally, the power 7^ is 

7„ = (5 - 3a)/(3 - a) , (51) 

with a being, we recall, the infoentropic index in Renyi statistics, which is limited to the 
interval 1 < a < 5/3. 

"Anomalous" diffusion is also called-for the explanations of the results of measurements 
in the case of experiments of cyclic voltammetry in microbatteries with fractal-like struc- 
tured electrodes As a consequence of the nowadays large interest associated to the 
development of microbatteries, the study of growth, annealing, and surface morphology of 
thin-films depositions used in cathodes, has acquired particular relevance 0]. These kind 
of systems are characterized by microroughnessed surfaces in a geometrically constrained 
region (nanometer scale), and then fractal characteristics can be expected to greatly in- 
fiuence the physical properties 5^. In cyclic voltammetry the resulting current between 
electrodes is determined by the charges arriving at the fractal electrode in a process of hy- 
drodynamic motion. Characterizing such motion as a Fickian diffusion fails to explain the 
experimental data, what can be done using a postulated law of "anomalous" diffusion 
The use of the latter is, as noticed before, a fitting procedure in a situation that requires 
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a higher-order thermo-hydrodynamics (motion of the fluid of charges in the nano- and sub- 
nanometric spaces around the surface of the fractal-hke structured electrode). The use of a 
zeroth-order thermo-hydrodynamics - the one leading to Fickian diffusion equations in a BG 
thermo-statistics - fails to satisfy Fisher's criterion of sufficiency, and if one persists in using 
a zeroth-order hydrodynamics, it needs be handled in an unconventional approach. Such 



31 



study has been performed resorting to Renyi statistics with details given in Refs. 
leading to Eq. (j^Uj) . Here we only noticed that, keeping fixed all physical characteristics 
in the experiment, if the morphology of the fractal electrode surface is modified, it can be 
established a relation between Renyi infoentropic index and the average fractal dimension 
df in the form 

a = {Adf - 1)1 {2df - 3) . (52) 

We recall that values of a are permitted in the interval 1 < a < 5/3 (otherwise the theory 
gives rise to singularities), now a clear physical restriction once then 2 < rfj < 3 as it should. 

C. Ideal Gas in Finite Box 

Let us now consider a simple but quite clarifying case, namely an ideal gas in a finite box 
(details in Ref. 0])- According to the exact result in R. K. Pathria textbook, jl^ and see 
also in the case of a finite but large box there follows for the energy per particle in the 
semiclassical limit the expression 



E 3, ^ 

— ~ -kgT 

N 2 



1 A\t 

1 + 



12 V ■ (^^' 

where V is the volume and A the area of the box, and At = ^/ VmkBT is the mean thermal 
de Broglie wavelength of the particles in the grand-canonical ensemble at temperature T. 
Evidently, in the thermodynamic limit it is recovered the result of energy equipartition. On 
the other hand, calculating in the thermodynamic limit but introducing the equivalent of 
the grand-canonical distribution in Renyi approach, in order to account for the insufficiency 
of description, we obtain that 



E 3^ 1 
N 2 



l + ^(a-l) 



(54) 
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where Fha is the Lagrange multipher associated to the energy and then, equating Eqs. 
and one finds that 



In the thermodynamic hmit A/V goes to zero and F^^ goes to ksT and then a goes to 1 as 
it should. Moreover, from Eq. ()55|) we can clearly see that the infoentropic index a depends 
on the system dynamics, its geometry and size, and the thermodynamic state. 

D. Black Body Radiation in Insufficient Description 

Finally, we briefiy mention the case of black-body radiation in insufficient description. 
We consider the gas of photons in the presence of an uniform fiux of energy (generated, 
for example, as a result of the presence of different temperatures at both ends of the con- 
tainer). We look for the calculation of the energy what is done, on the one side, using the 
conventional approach in a description that includes the energy and the energy flux as ba- 
sic variables (that is, in a flrst-order t her mo-hydro dynamics; see discussion on higher-order 
thermo-hydrodynamics earlier in this Section) and, on the other side, using an incomplete 
description including only the energy (zeroth-order thermo-hydrodynamics) but dealt with 
using Renyi statistics trying to compensate for the insufficiency of description. Details are 
given in Ref. jsil, and here we only notice that equating the values of the energy in both 
descriptions, it follows that the infoentropic index is given by 



where / is the energy flux and Jq = (8/15y^)a^c^T^ with a being Stefan- Boltzmann constant. 
It can be noticed that Jq is a kind of a flux of energy composed of the energy density of 
the radiation, aT^ /V , traversing at the speed of light, while I is roughly the density of 
energy traversing at a speed determined by the gradient of temperature, and then ///q is 
very small, and then the infoentropic index is practically 1 meaning that the insufficiency 
in the characterization can be ignored. 




(55) 



a = 1 - (///o)^ 



(56) 
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IV. CONCLUDING REMARKS 



However the enormous success and large application of Shannon- Jaynes approach to 
derive the probability distributions of the ensemble formalism in the Laplace-Maxwell- 
Boltzmann-Gibbs statistical foundations of physics, as it has been noticed, some cases are 
difficult to be properly handled within such formulation, as a result of existing some kind 
of fuzziness in data or information, that is, the presence of conditions of insufficiency in the 
characterization of the (microscopic, macroscopic or mesoscopic) state of the system. Such 
difficulty with the proper characterization of the system in the problem in hands can be 
somehow compensated, as shown, with the introduction of peculiar parameter-dependent 
alternative informational-entropies (see Table I), allowing for the construction of a vast 
number of auxiliary heterotypical statistics. In that way one recovers the ability to make 
improved predictions on the properties of the system, and he capable to get a picture of the 
physical processes involved allowing, as for example illustrated in Section III, to evaluate 
the influence of peculiar characteristics of the sample (like, e.g., its fractal structure) on its 
physico-chemical behavior. 

Returning to the question of the possible failure to satisfy Fisher's criteria of efficiency 
and sufficiency, in the case of particularly dealing with systems with some kind of fractal-like 
structure the use of Shannon- Jaynes infoentropy would require to introduce as information 
the highly correlated conditions that are in that case present. Two examples in condensed 
matter physics have been described in the previous Section in which fractality enters via 
the non-smooth topography of the boundary surfaces which have large influence on phe- 
nomena occurring in constrained geometries (nanometer scales in the active region of the 
sample). Hence the most general and complete Boltzmann-Gibbs formalism in Shannon- 
Jaynes approach becomes hampered out and is difficult to handle, and then, as shown, use 
of other types of informational-entropies leads to the derivation of heterotypical probability 
distributions on the basis of the constrained maximization of unconventional informational- 
statistical entropies (quantity of uncertainty of information), to be accompanied, as noticed 
in the main text, with the use of the so-called escort probabilities, allowing for an analysis 
and understanding of the observed physical behavior of the system. 

Summarizing, Unconventional Statistical Mechanics consists of two steps: 1. The choice 
of a deemed appropriate structural informational-entropy for generating the heterotypical 
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statistical operator, and 2. The use of a escort probability in terms of the heterotypical 
distribution of item 1. 

As shown in the main text, and illustrated in Section III, the escort probability introduces 
corrections to an inefficient description by including correlations and higher-order variances 
of the observables involved. On the other hand, the heterotypical distribution introduces cor- 
rections to the insufficient description by modifying the statistical weight of the dynamical 
states of the conventional approach involved in the situation under consideration. More- 
over, we have considered a particular case, namely the statistics as derived from the use of 
Renyi informational entropy, centering the attention on the derivation of an Unconventional 
Statistical Mechanics appropriate for dealing with far-removed-from-equilibrium systems. 
Moreover, we have reported the calculation, in such conditions, of the distribution functions 
of single fermions and bosons, the counterpart in these unconventional statistics of the usual 
Fermi-Dirac and Bose-Einstein distributions: These distributions are illustrated in Figs. 
and El 

In conclusion, we may say that USM appears as a valuable approach, in which the 
introduction of informational-entropic-indexes-dependent informational-entropies leads to 
a particularly convenient and sophisticated method for studying certain classes of physical 
systems for which the criteria of efficiency and sufficiency in its characterization cannot be 
properly satisfied. 

Finally, it is relevant to emphasize again the fact that the infoentropic index(es) is(are) 
dependent on the dynamics involved, the system's geometry and dimensions, boundary 
conditions, its macroscopic-thermodynamic state (in equilibrium, or out of it when becomes 
a function of time), and the experimental protocol. 

Also, we recall the fundamental point, stressed in the Introduction, that these informa- 
tional entropies (informational measures) are not at all to be interpreted as the physico- 
thermal entropy of the system. 
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FIG. 1: The distribution of Eq. (jSHJ for a kinetic temperature [cf. Eq. (42)] of 300K and values 
of Renyi's infoentropic- index a smaller than 1. 
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FIG. 2: The distribution of Eq. (jSHJ for a kinetic temperature [cf. Eq. (42)] of 300K and values 
of Renyi's infoentropic- index a larger than 1. 
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FIG. 3: The high-energy side of the photoluminescence spectra for the quantum-well widths indi- 
cated in the upper right inset. The continuous lines are the best fittings to experimental curves 
using Eq. (|^ . 
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FIG. 4: Dependence of the infoentropic index a with the quantum-well width following the ap- 
proximate empirical law of Eq. (|48() . 



33 



